Role of spatial coherence in polarization tomography 
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We analyze an experimental setup in which a quasi-monochromatic spatially coherent beam of 
light is used to probe a paraxial optical scatterer. We discuss the effect of the spatial coherence of 
the probe beam on the Mueller matrix representing the scatterer. We show that according to the 
degree of spatial coherence of the beam, the same scattering system can be represented by different 
Mueller matrices. This result should serve as a warning for experimentalists. 
OCIS codes: 030.1640, 260.5430, 290.0290. 
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A great deal of literature exists about the phenomeno- 
logical treatment of polarization of light by matrix meth- 
ods 0. Among these methods, the Mueller-Stokes 
one is perhaps the most utilized for the description of 
the interaction between a quasi-monochromatic beam of 
light, hereafter denoted as the probe, and a polarization- 
affecting medium, namely the scattering system. In the 
Mueller-Stokes formalism the beam of light is represented 
by four real numbers, namely the Stokes parameters, and 
the scattering system is represented by a 4 x 4 real ma- 
trix, the Mueller matrix. Recently [2j, considerable ef- 
forts have been devoted to the study of the connection 
between the algebraic properties of Mueller matrices and 
the corresponding physical properties of the systems they 
represent. However, in most contributions [ij, |2| empha- 
sis is given to the mathematical aspects of the problem 
rather than to the physical ones. 

On the contrary, in the present Letter we study how 
the spatial coherence of the probe beam affects the phys- 
ical properties of the Mueller matrix representing a given 
scattering system. Using Wolf's unified theory of coher- 
ence and polarization of random electromagnetic beams 
, we demonstrate that the same scattering system can 
be represented by different Mueller matrices, the differ- 
ence depending on the degree of spatial coherence of the 
probe beam. Specifically, we show that an optical scat- 
terer may behave either as a non-depolarizing or a de- 
polarizing system according to whether the probe beam 
was completely spatially coherent or completely spatially 
incoherent. 

A typical experimental setup for polarization tomog- 
raphy consists of five elementary units: The source S, 
the polarization-preparer V, the scattering system M 
the polarization-analyzer A, and the detector V (see Fig. 
1). The source S emits a quasi-monochromatic beam of 
light at frequency u>, which is prepared by V in a well 
defined pure polarization state. This beam propagates 
close to the z-axis through a scattering system charac- 
terized by its spectral transmission matrix Gy(ri, ro, w), 
(i = x,y, j = x,y), where r = (po,z ), r 1 = (pi,zi), 
and pa = (xa,Ua), (A = 0,1), are the transverse posi- 
tion vectors on the input (z = Zq) and output (z = z\) 



planes, respectively. Let Ej(ro,oj), (j = x,y) be a spec- 
tral Cartesian components of the electric field at the 
point Qo of coordinates ro = (po, zq), in the input plane. 
The Cartesian components of the electric field at the 
point P\ of coordinates ri = (pi, z\) in the output plane 
are then given by: 



E. l (r 1 ,uj) = / d 2 po G ij {ri,r ,uj)E j 



(ro,w), (1) 



(i = x,y; j = x,y), where the integration extends over 
the transverse input-plane coordinates: d 2 po = dzodyoi 
and summation on repeated indices is understood. Be- 
hind the scatterer, the polarization-analyzer A can select 
an arbitrary field component E(yi,oj) = a x E x (r\,uS) + 
CLyEy (ri , oj) , where a Xl a y are complex amplitudes de- 
termined by the analyzer setup. Finally, the detector 
V detects the spectral density S(t\,uj) (proportional to 
the averaged electric field density) of the beam at the 
point P\. S{t\,uj) = (E(ri,uj)E*(ri,uj)), where angular 
brackets denote the average over an ensemble of real- 
izations of the scattered electric field 0- This expres- 
sion for S(ri,u>) is a particular case of the more general 
formula S(r,cu) = Tr{W(r,w)} 6], where Tr{W(r,w)} 
denotes the trace of the 2x2 spectral density matrix 
W(r, lo) = W(r, r, cu) evaluated at the point P (specified 
by the position vector r), and \N(ta,tb,u>) is the cross- 
spectral density matrix of the beam at the pair of points 
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FIG. 1: Polarization tomography setup described in the text. 
The input (z = Zo) and the output [z = z\) planes are indi- 
cated by vertical lines. 
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?A = (Pa, Zq), r B = (pa, zo) in the plane z = z : 

W i:i (r A ,r B ,u) ee (E t (r A ,uu)E*(v B ,uu)), (2) 

where (i = x, y; j = x,y) ■ 

We now consider the effects of the propagation through 
the scattering system on the spectral density matrix 
W(r,w) of the beam. From Eqs. (QJ and J5J it read- 
ily follows that 

W(r x ,w) = | d 2 p d 2 p '{ G ( r i> r o^) (3) 
xW(4ri»G(r;',r llW )}, 

where summation on repeated indices is understood and 
Wki(r' ,r%,u) ee (E k (r' ,ui)Ef(4,u))) are the elements 
of the cross-spectral density matrix W(r ,r ',cj) of the 
beam at all pairs of points p' , p ' in the input plane, and 
r o = (Poi z o): r o = (Po' z o) respectively. Furthermore, 
Gij(rQ,Ti,u) = G*j(r!, r ', lu). In order to obtain a for- 
mula involving Stokes parameters and Mueller matrices, 
we multiply both sides of Eq. © by the normalized Pauli 
matrix ai a -\ (a = 0, 1, 2, 3) |8( and trace, obtaining 



S a (r 1 ,uj) 

= f d 2 P ' d 2 ^^,^,4,^(^,4,0,), 



(4) 



(f3 = 0,1,2,3), where summation on repeated indices 
is understood and we have defined the spectral density 
Stokes parameters of the beam in the output plane 
as: S a (ri,cj) ee Tr{<7( Q )W(ri, w)}. Moreover, we have 
introduced the cross-spectral density Mueller matrix 



ee Tr{cr (ct) G(ri, r' , w)cr (/3) G t (ri, r ', u>)}, 



(5) 



and the cross-spectral density Stokes parameters of the 
input beam 



J (r' o ,r^)=Tr{a w W(r' o ,4,cj)}, 



(6) 



which reduces to the input spectral density Stokes pa- 
rameters Sp(r' ,uj) in the limit r ' — > r . 

Eq. J3J shows that unlike the Stokes parameters {S^} 
in the elementary theory of partial polarization, the spec- 
tral density Stokes parameters {S a (ri, u>)} at any point 
Pi in the output plane are not just a linear combi- 
nation of the input spectral density Stokes parameters 
{S a (ro, uj)} but they are expressed in terms of the cross- 
spectral density Stokes parameters { Jp(r' , r ', uj)} at all 
pairs of points ppiPo m the input plane. This fact ob- 
viously represents the effect of the coherence of the in- 
put beam. Moreover, since for r ^ r ', G^(ri, r ', lu) ^ 
[G(ri,r ,o;)]^, the cross-spectral density Mueller matrix 
is not a Mueller- Jones matrix Q. 

So far we have been concerned with a generic spatially 
coherent light source, without specifying its degree of co- 
herence. We shall now consider with some detail the 



two opposite limit cases of completely coherent and com- 
pletely incoherent light. For both cases we consider a uni- 
formly polarized input beam specified by the coordinate- 
and frequency-independent two-dimensional unit vector 
e = (e x ,e y ), such that Ei(r Q ,uj) = eiE(r ,u), (i = x,y), 
where E(tq, lu) is a scalar function of the point Qq in the 
input plane. Using this assumption, one readily obtains 
the following expression for the cross-spectral density ma- 
trix of the input beam: 



W k i(r' ,r'o,u)) = £ k iw(r' ,r'o,uj), 



(7) 



where £ M = e k e\, and w(r' Q , r ', lu) ee (E(t' q ,lu)E* (y'^lu)) 
is the scalar cross-spectral density function which charac- 
terize the second-order coherence properties of the input 
beam. If we substitute Eq. into Eq. @ and use Eqs. 
(I5I6H we obtain 



(8) 



S Q (r l7 Lu) = J dVod 2 p Xr ,r >) 

XM a /3(Tl,T' o ,4,Uj)S , 



where Sp = Tr{atp\£}, (f3 — 0,1,2,3), are the input 
Stokes parameters which are independent from both the 
input-plane coordinates ro and the frequency lu. 
In the case of a completely coherent source 



w(r ,r ',w) =u(r' ,co)u*(4,uj), 



(9) 



where u(ro,w) represents the complex disturbance of the 
field in the input plane [13 . From Eqs. © and © read- 
ily follows that S a (ri , lu) = M^p (ri, u)Sp where the "co- 
herent" Mueller matrix M^ c '(ri,u>) at the output point 
Pi is defined as 

M%?(t u u) ee U{ ff(a) K(r ll uV (w K , (ri ) «)} 1 (10) 

where K(ri,u>) ee J d 2 po w(r , w)G(ri, r , uj) represents 
the tensor- valued complex disturbance of the field in the 
output plane. Eq. (|1C)|> shows that in the ideal case 
of a completely coherent and uniformly polarized probe 
beam, the Mueller matrix M^(v\, lu) is a Mueller- Jones 
matrix Q representing a non- depolarizing scattering sys- 
tem. 

In the case of a completely incoherent source 

w(r' ,4,cu) = w(r' ,uj)5W(p' -pZ), (11) 

where w(r' ,uj) > is the spectral density function of 
the input beam. If we substitute Eq. 111(1 in Eq. (jHJ) 
we obtain S a (ri,ui) = (ri, uu)Sp where the "inco- 
herent" Mueller matrix M ^ (ri , lu) at the output point 
Pi is defined as 

Mg(n,u;)= / ' d 2 p' w(v' ,uj)M a p(r 1 ,r' ,uj), (12) 

where the spectral density Mueller matrix 
M a p(Y\, r , uj) ee M a p(rx, r , r , uj) is now a bona 
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fide Mueller-Jones matrix. Eq. (|12f> expresses the 
incoherent Mueller matrix M V> (ri , ui) as a linear com- 
bination, with non-negative coefficients w(r' Q ,u>), of 
the Mueller-Jones matrices A4 a g(ri,r' Q ,u)). Therefore 
M ^ (ri , uj) is formally equivalent to a Mueller matrix 
representing a depolarizing scattering system 0. More- 
over, from Eq. 1(12(1 readily follows that the Hcrmitian 
matrix H 11] associated with M^(ri, uj) is positive 
scmidcfinite, namely M^(ri,w) satisfy the so called 
Jones criterion and can be referred to as a "physical" 
Mueller matrix |2|. 

Eqs. (|lt)(l and 1)12(1 are the key results of this Let- 
ter: they give the Mueller matrices M (ri , uj) and 
M^(ri, uj) describing the same scattering system probed 
by either a completely coherent or a completely incoher- 
ent light beam, respectively. Surprisingly, we have found 
that in the first case M^^ri, uj) is a Mueller- Jones ma- 
trix which represents a non-depolarizing optical system, 
while in the second case M^(ri, uj) is a Jones-derived (or 
"physical" ) Mueller matrix that represents a depolarizing 
optical system. 

The more general (and realistic) case of partially co- 
herent quasi-monochromatic light probe is illustrated 
by Eqs. 141618(1 . In particular Eq. © shows that 
when a partially coherent probe beam is used, the 
cross-spectral density Mueller matrix may not be inter- 
preted as a Mueller- Jones matrix. However, since in 
general the cross-spectral density function can be ex- 
panded as a superposition of coherent fields {u n (po, w)} 
as w(r' ,r%,u) = ^ n -f n (z ,u)u n (p' ,ui)ul(p^,uj) 0, 
where 7„(z ,w) > 0, it is easy to see that the "partially- 
coherent" Mueller matrix M^ p '(ti,u>) 

M<$(t u w) = J d 2 p' d 2 rtw(r> Q ,r'>,Lo) ^ 
xX Q(3 (n,ro,ro,cj), 

can be also decomposed as a superposition 
[with non-negative coefficients r y n (zQ,uj)] of "co- 
herent" Mueller- Jones matrices M^(ri,o;) = 
Tr{a( Q )K„(ri,w)cr (/3) K,J l (ri,w)} where K n (ri,a>) = 
J d 2 p u n (r , w)G(i"i, r , uj). Therefore the matrix 
M^ P \r i 7 uj) represents, in general, a depolarizing optical 
system. It is worthwhile to note that each matrix 
M^(ri,o;) depends in a nontrivial way on the field- 
amplitudes u n (ro,uj) of the probe beam. Therefore, 
when using spatially coherent probe beams, the linearity 
required by the Mueller-Stokes formalism can be ensured 
only with respect to the polarization degrees of freedom 
of the beam (viz, the matrix indices), but not with 
respect to the spatial degrees of freedom. 

In conclusion, we have shown that, contrary to com- 
mon belief, the Mueller matrix representing an optical 
system in a scattering experiment, is not determined by 
the system only but strongly depends on the spatial co- 



herence of the probe beam. This fact poses serious limi- 
tations to the range of applicability of the Mueller-Stokes 
formalism. Consider, for instance, the following case: a 
probe beam propagates through a medium suffering mul- 
tiple scattering. Since the degree of spatial coherence of 
the beam changes during propagation 0, the last part 
of the medium will be probed by a beam qualitatively 
different from the one which probed the first part of the 
medium. Therefore the Mueller matrix representing the 
medium will not only be determined by the polarization 
properties of the medium itself, but also by the way the 
medium affects the coherence properties of the probe 
beam. With respect to this problem, the use of com- 
pletely spatially incoherent probe beams for polarization 
tomography seems highly preferable. 
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